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Abstract 

We present an exact calculation of S and P wave QQnn states using different standard non- 
relativistic quark-quark potentials. We explore in detail the charm and bottom sectors looking 
for bound states that could be measured within existing facilities. Against the proliferation of 
four-quark states sometimes predicted in the literature, we found a small number of candidates 
to be stable. We analyze their properties in a trial to distinguish between compact and molecular 
states. Possible decay modes are discussed. 

PACS numbers: 12.39.Jh,14.40.Lb,21.45.+v,31.15.Ja 
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I. INTRODUCTION 



The existence of stable QQnn states has been the topic for discussion since the early 
80's pj. These states are of particular interest since they are manifestly exotic, i.e., heavy 
flavor quantum number ±2 with baryon number equal 0. If they would lie below the thresh- 
old for dissociation into two ordinary hadrons they would be narrow and should show up 
clearly in the experimental spectrum. There are already estimates of the production rates 
indicating they could be produced and detected at present (and future) experimental fa- 
cilities [2J. Four-quark states seem to be necessary to tame the bewildering landscape of 
new meson states. Non-exotic four-quark states embedded in the meson spectra have been 
proposed as a thoughtful explanation of the proliferation and peculiar properties of light 
scalar-isoscalar or open-charm mesons [3j. There are recent indications of the possible ex- 
istence of molecular four-quark states, suggested from the observation of the Z + (4430) 
A step forward scrutinizing the structure of low-energy hadrons would be to demonstrate 
the existence of stable manifestly exotic multiquark states, the QQnn system being an ideal 
candidate. 

To illustrate the theoretical landscape, we present in Table II a summary of different ap- 
proaches to the manifestly exotic four-quark spectroscopy (1, Ej, 0, H, @, Eo, [0, E, 13, 0, 



151 . Il6l . 1171 ] . Exotic multiquarks were examined in Ref. [5J] solving the four-body problem 
by three different variational methods with a nonrelativistic potential considering explicitly 
virtual meson-meson components in the wave function. In Ref. [6j four-quark states were 
studied using a variational approach with trial wave functions whose interaction energies 
are approximately given by known hadron masses. Ref. Q used a potential derived from 
the MIT bag model in the Born-Oppenheimer approximation. The calculations were done 
by means of the Green's function Monte Carlo method, bn — bn molecules loosely bound b 
the one-pion exchange were obtained in Ref. using chiral perturbation theory. Ref. 
analyzed L = four-quark systems with the Bhaduri potential and a variational method 
in a harmonic oscillator basis up to N = 8 quanta. Ref. [loj discussed the stability of 
multiquark systems using different parametrizations of a Goldstone boson exchange model 



and a variational formalism with gaussian trial wave functions. In Ref. [llj the Bhaduri 
potential was reexamined by means of a variational method that allows nonzero internal 
orbital angular momentum in the subsystems of quarks and antiquarks. The existence of a 
shallow tetraquark state ccud was discussed in Ref. [l2[ using semiempirical mass relations. 
Ref. [13} designed a powerful method, similar to the stochastic variational approach (l8| . 
accommodating two free-meson asymptotic states. It was applied to the Bhaduri potential, 
fixing the results of Refs. [j^ 11]. Ref. HI analyzed multiquark states with a variational 



formalism using gaussian trial wave functions with only quadratic terms in the Jacobi co- 
ordinates. QCD sum rules were used in Re f. |l5l| to study the possible existence of an axial 
diquark-antidiquark bound state. In Ref. [16| the ground state of tetraquarks were evalu- 
ated assuming a diquark-antidiquark structure, reducing the relativistic four-body problem 



to the solution of two relativistic two-body problems. Ref. [17( discussed the possible exis- 
tence of four-quark bound states within the framework of the chiral SU(3) quark model by 
means of a variational approach using gaussian trial wave functions. 

As seen in Table HJ there is a remarkable agreement on the existence of a / = 0, J p = l + 
bbnn bound state and also, although not so neatly, on the existence of a ccnn one (along this 
work n stands for a light u or d quark). However, among the different theoretical approaches 
only a few had paid attention to other quantum numbers, trying to elucidate if a proliferation 
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of four-quark states is predicted. Hence, in this work our purpose will be twofold. On the 
one hand, we shall try to shed some light on this topic by making a detailed analysis of the 
QQnn spectra discussing each set of quantum numbers. On the other hand, the Achilles' 
heel of almost all the approaches described above is the lack of an exact numerical method 
to solve the four-body problem. Most of these works rely on variational calculations with 
different types of trial wave functions or on semiempirical mass relations. The importance 
of this problem requires numerical methods able to provide exact solutions with controlled 
numerical uncertainties. 

To this end we shall use a new approach based on the hyperspherical formalism recently 



developed to solve exactly the four-quark problem |19l . I20j ]. The idea is to perform an 
expansion of the trial wave function in terms of the hyperspherical harmonics (HH) functions. 
This allows to generalize the simplicity of the spherical harmonic expansion for the angular 
functions of a single particle motion to a system of particles by introducing a global length 
p, the hyperradius, and a set of angles, Q. For the HH expansion to be practical, the 
evaluation of the potential energy matrix elements must be feasible. The main difficulty of 
this method is to construct HH functions with proper permutational symmetry for a system 
of identical particles. This is a difficult problem that may be overcome by means of the HH 
formalism based on the symmetrization of the N— body wave function with respect to the 
symmetric group using the Barnea and Novoselsky algorithm j2l[. For systems containing 
only two pairs of identical particles this problem is greatly simplified making the appropriate 
choice of Jacobi coordinates. In Ref. [20] we have developed the HH formalism for the QQnn 
system, and we will use it here with the appropriate modifications to study the exotic QQnn 
four-quark spectra. 

The manuscript is organized as follows. In Sec. [IT] we briefly revise the HH formalism for 
the QQnn system and present the quark models used. In Sec. II III we introduce observables 
that may allow to distinguish between unbound and compact or molecular four-quark bound 
states. In Sec. IIVI the results and the analysis of the ccnn and hhnn spectroscopy are 
presented. In Sec. |V]the possibility of measuring the predicted bound states within current 
experimental facilities is discussed. Finally, we summarize in Sec. I VI I our conclusions. 



II. TECHNICAL DETAILS 

Within the HH expansion, the four-quark wave function can be written as a sum of outer 
products of color, isospin, spin and configuration terms 

\4>cisr) = |Color)|Isospin) [|Spin) ® \R)] JM , (1) 

such that the four-quark state is a color singlet with well defined parity, isospin and total 
angular momentum. In the following we shall assume that particles 1 and 2 are the Q-quarks 
and particles 3 and 4 are the n-quarks. Thus, particles 1 and 2 are identical, and so are 3 
and 4. Consequently, the Pauli principle leads to the following conditions, 

P\2\<Pcisr) = Pu\<Pcisr) = -\<Pcisr) , (2) 

Pij being the permutation operator of particles % and j. 

Coupling the color states of two quarks (antiquarks) can yield two possible represen- 
tations, the symmetric 6-dimensional, 6 (6), and the antisymmetric 3-dimensional, 3 (3). 
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Coupling the color states of the quark pair with that of the antiquark pair must yield a color 
singlet. Thus, there are only two possible color states for a QQqq system j22| . 

|Color> = {|3 12 3 34 >, |6i 2 6 34 )}. (3) 

These states have well defined symmetry under permutations, Eq. ([2]). The spin states with 
such symmetry can be obtained in the following way, 

|Spin) = |(0i, s 2 )5'i2, (s^Si)S M )S) = \(S 12 Su)S) . (4) 

The same holds for the isospin, |Isospin) = | (z 3 , z 4 )/ 3 4) , which applies only to the n-quarks, 
thus I = J34. 

We use the HH expansion to describe the spatial part of the wave function. We choose 
for convenience the if-type Jacobi coordinates, 



ri! = A*i, 2 (r* 2 - ri) , 

'm 3 r 3 + m A r A min + m 2 r 2 \ 

m 34 77112 / 

% = /Mr^ - r s) , (5) 



f] 2 = A*12,34 



were my = mi+rrij, p,^ = ^rriirrij / , and 777.1234 = 7771+777.2+7773+7774. Using these vectors, 
it is easy to obtain basis functions that have well defined symmetry under permutations of 
the pairs (12) and (34). In the HH formalism the three Jacobi vectors are transformed into 
a single length variable, p = ^/r/f + 77! + 773 , and 8-angular variables, f2, that represent the 
location on the 8-dimensional sphere. The spatial basis states are given by 

( P n\R) = u n ( P )y [K] (n) , (6) 

were y\K\ are the HH functions, and [K] = {KK 12 LM L L 12 l 3 d 2 (ii\. The quantum number K 
is the grand angular momentum, LM^ are the usual orbital angular momentum quantum 
numbers, and i{ is the angular momentum associated with the Jacobi vector 77^. The quan- 
tum numbers K\ 2 , L\ 2 correspond to the intermediate coupling of rj 1 and f] 2 - The Laguerre 
functions are used as the hyper-radial basis functions U n (p). 

The Pauli principle, Eq. (j2j), leads to the following restrictions on the allowed combina- 
tions of basis states: 

(i) (-l) s i2+ £ i = +1, (_i)S 3 4+/+^ = _i for t h e |6 12 6 34 ) color state, 

(ii) (-1)^2+^ = -1, (_i)S 3 4+/+£ 3 = + l for the |3i 2 3 34 ) state. 

Assuming non-relativistic quantum mechanics we solve the four-body Schrodinger equa- 
tion using the basis states described above. The grand angular momentum K is the main 
quantum number in our expansion and the calculation is truncated at some K value. 

For our study we will use two standard quark potential models. The constituent quark 
cluster (CQC) model was proposed in the early 90 's in an attempt to obtain a simultaneous 



description of the nucleon-nucleon interaction and the baryon spectra [23j. Later on it was 
generalized to all flavor sectors giving a reasonable description of the meson [24] and baryon 
spectra (25|. The model contains Goldstone boson exchange between quarks, a one-gluon- 
exchange (OGE) potential, and a screened confined interaction as dictated by unquenched 



lattice calculations [26|. The model parameters have been taken from Ref. [24j with the 
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FIG. 1: H-type Jacobi vectors. 




exception of the OGE regularization parameter, see Ref. 20| for details. In the following we 
shall denote this parametrization as CQC and the standard parametrization of Ref. 24j will 
be referred to as CQCis- Explicit expressions of the interacting potentials and a detailed 



discussion of the model can be found in Ref. [24(. The Bhaduri-Cohler-Nogami (BCN) model 



was proposed in the early 80's in an attempt to obtain a unified description of meson and 
baryon spectroscopy (2?|]. It was later on applied to study the baryon spectra [28| and four- 
quark systems The model retains the most important terms of the one-gluon exchange 
(OGE) interaction, namely coulomb and spin-spin, and a linear confining potential. The 
parameters are taken from Ref. [9|. 

A summary of the energies obtained with both models, CQC and BCN, for selected 
meson states that may appear in the thresholds of the studied four-quark systems are given 
in Table [IIJ together with the corresponding experimental energies. 



III. BOUND STATES, COMPACT AND MOLECULAR 
A. Threshold determination 

As thoroughly discussed in Ref. 20|, in order to discriminate between four-quark bound 
states and simple pieces of the meson-meson continuum, one has to carefully determine the 
two-meson states that constitute the thresholds for each set of quantum numbers. Dealing 
with strongly interacting particles, the two-meson states should have well defined total an- 
gular momentum (J), parity (P), and a properly symmetrized wave function if two identical 
mesons are considered (spin-statistics theorem). When noncentral forces are not taken into 
account, orbital angular momentum (L) and total spin (S) are also good quantum numbers. 
We give in Tables 1111} HV] |Vj IVIj and IVHI the lowest threshold in both cases, that we will 
refer to as coupled (CO) and uncoupled (UN) schemes respectively, together with the final 
state relative orbital angular momentum of the decay products. We would like to emphasize 
that although we use central forces in our calculation the CO scheme is the relevant one for 
observations, since a small non-central component in the potential is enough to produce a 
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sizeable effect on the width of a state. In Table II I II we summarize the thresholds obtained 
using the experimental energies given in Ref. 29]. In Tables HVl and IVl we quote the thresh- 



olds obtained with the CQC model for the charm and bottom sectors respectively, and in 
Tables EH and EH] those for the BCN model. 

An important property of the QQnn system, that is crucial for the possible existence of 
bound states, is the fact that only one physical threshold (Qn)(Qn) is allowed. Consequently, 
particular modifications of the four-quark interaction, for instance a strong color-dependent 
attraction in the QQ pair, would not be translated into any asymptotically free two-meson 
state. As discussed in Ref. |31|, this is not a general property in the four-quark spectroscopy, 



as the QQnn four-quark state has two allowed physical thresholds: (QQ)(nn) and (Qn)(nQ). 



B. Figures of merit 

The relevant quantity for analyzing the stability of any four-quark state is Ae, the energy 
difference between the mass of the four-quark system and that of the lowest two-meson 
threshold, 

A E = E 4q -E(M 1 ,M 2 ), (7) 

where E± q stands for the four-quark energy and E{M\, M 2 ) for the energy of the two-meson 
threshold. Thus, Ae < indicates that all fall-apart decays are forbidden, and therefore one 
has a proper bound state. A E > will indicate that the four-quark solution corresponds to 
an unbound threshold (two free mesons). 

One of the main difficulties in studying four-quark states, already discussed in Refs. [ljj 



201 ]. is the slow convergence of unbound solutions. We show in Fig. [2] the evolution of 
E± q /E(M\, M 2 ) as a function of K for three different states, two of them bound (solid and 
dash-dotted lines) and one unbound (dashed line). Although all of them converge for large 
enough values of K, the correct description of the two-meson threshold for unbound states is 
slow and time consuming. A helpful tool to minimize this problem was proposed in Ref. 19l 
through an extrapolation of the four-quark energy using the expression 

E 4q (K) = E 4q (K = 00) + -^, (8) 

where E iq (K = 00), a and b are fitted parameters. When this extrapolation is used for 
unbound states, one can observe how the four-quark energies reproduce the thresholds to 
within a few MeV. 

A second important quantity to characterize the possible existence of a bound state is 
the evolution of the root mean square radius as a function of K. While for unbound states, 
the components of the four-quark state will tend to be far away when increasing K, for a 
bound state the radius should stabilize when increasing K. In order to compare four-quark 
against two free-meson states we define the root mean square radius for four (two) quark 
systems 

'T£ lmi({ri- Rf) 

e£? m i 



RMS 4q{2(?) = [ ^ =1 ^ ^ ] • (9) 



and its corresponding ratio 



= RMS 4l? 

R RMS Mx + RMS M2 ' 1 ; 
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FIG. 2: E 4q /E(M 1 ,M 2 ) as a function of K for the J P (L,S,I) = l+(0, 1,0) (solid line), 0+ (0,0,0) 
(dashed line), and 1 _ (1,0, 1) (dash-dotted line) with the CQC model. 




where RMSmi + RMSm 2 stands for the sum of the radii of the mesons corresponding to the 
lowest threshold. 



C. Compact vs molecular states 

Besides trying to unravel the possible existence of bound QQnn states one should try 
to understand whether it is possible to differentiate between compact and molecular states. 
A molecular state may be understood as a four-quark state containing a single physical 
two-meson component, i.e., a unique singlet-singlet component in the color wave function 
with well-defined spin and isospin quantum numbers. One could expect these states not 
being deeply bound and therefore having a size of the order of the two-meson system, i.e., 
Ar ~ 1. Opposite to that, a compact state may be characterized by its involved structure on 
the color space, its wave function containing different singlet-singlet components with non 
negligible probabilities. One would expect such states would be smaller than typical two- 
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meson systems, i.e., A# < 1. Let us notice that while Ar > 1 but finite would correspond 
to a meson-meson molecule A# K .Zl°? qq would represent an unbound threshold. 

One may illustrate the situation described above by the deuteron and the if-dibaryon 
examples. Let us try to draw the analogy between these states and the QQnn system. The 
deuteron has a small binding energy of —2.225 MeV and a ratio between its root mean square 
charge radius (2.139 fm) and the one of two protons (1.75 fm) of 1.222 [32|], A# = 1.22 and 
A E = —2.225 MeV in our notation. Should the deuteron be considered as a pure baryon- 
baryon molecule? Although Ref. [33] emphasized the difficulties to identify pure hadron- 
hadron molecules close to thresholds, in the deuteron case it was long-ago justified [Hj]. The 
probability of physical two-baryon states other than nucleon-nucleon is meaningless (23| . 
Therefore it constitutes a clear example of a molecular state. The postulated if-dibaryon 
would however fit into the picture of compact states, its wave function presenting relevant 
components of different singlet- singlet physical channels: AA, AH, and EE at least [351 ] . 

This last aspect makes contact with the role played by hidden-color configurations, color 
singlets built by nonsinglet constituents. There are three different ways for coupling two 
quarks and two antiquarks in a colorless state, 

[(?i?2)(M4)] = {|3i 2 3 34 >, |6 12 6 34 >} 
[(gi&XMOj = {II13I24), |8i 3 8 24 )} 

[(M4)(<M3)] = {|ll4l23>,|8 14 8 23 >}. (11) 

Each coupling scheme allows to define a color basis where the four-body problem can be 
solved. The first basis, [(gi^fe^)], being the most suitable one to deal with the Pauli 
principle, is made entirely by hidden-color vectors. The other two are hybrid basis that 
contain both singlet-singlet (physical) and octect-octect (hidden-color) components. It 
is possible to prove from simple group theory arguments that once we have solved the 
four-body problem for a system composed of two identical quarks (QQ) and two identical 
antiquarks (fin), there is a minimum value for the octect-octect component probability of the 
wave function either in the [(Qin 3 )(Q2n^)] or the [{Qin4)(Q 2 n 3 )} couplings: -Pgg' 2 , Pgg' 23 £ 
[1/3, 2/3]. It can also be proved that for a four-quark threshold state -Pgg 3 ' 24 = Pgg 4 ' 23 = 4/9. 
Does this imply an important hidden-color component in all QQnn states? The answer is 
no. In Ref. [36] it was proved than any physical state can be expanded in terms of a basis 
constructed by direct product of mesonic and/or baryonic states not necessarily linearly 
independent. Therefore one can express any QQnn state in terms of the singlet-singlet 
component of the [(Qin3)(Q2?k)] and [(Q 1^4) (Q 2113)) basis. 

This discussion can be made more quantitative. Let us assume that {P, Q} and {P, Q} 
are the projectors associated to two orthonormal basis that are not orthogonal to each other, 
i.e., PP I <p) 7^ and PQ | 0) 7^ for an arbitrary state | 0). This would be the case of the 
two orthonormal basis: { | I13I24), I813824)} and { | I14I23) , I814823)}. Any arbitrary state can 
be written as 

I *) = P I + Q I (12) 
and the probability of the state associated to P or P will be given by [381 ] . 

?!*>([«]) 



* I PQ I m ) + ( * I QP I 



2(1 — cos 2 a) 

= 2(1 -La) [<* I PQ I *) + (* I QP I *)] ' < 13) 
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where P = | u) (u | and P — \u') (u 1 | and cosa = (u' \ u). For a molecular state either 
V m {[u\) or V^\[u'}) would be equal to zero while for a compact state both will be different 
from zero. 



IV. RESULTS 

A. Comparison to other numerical methods 

To illustrate the performance of the numerical procedure described in Sec. [Ill it is conve- 
nient to compare with other numerical methods to understand its capability and advantages, 
if any. As outlined in Sec. [Q in the past decades there have been several attempts to study 
multiquarks containing explicit charm or bottom flavors. Among them we shall analyze the 
calculation of the (L, S, I) = (0, 1, 0) ccnn state of Refs. (§, [l3| using the BCN potential and 
the results obtained in Ref. [14J with the CQC model. 

We present in Table IVIIII results for different L = spin-isospin ccnn states within the 
CQC model. We quote in the first column the results obtained with a variational calculation 



using gaussian trial wave functions with only quadratic terms in the Jacobi coordinates 14 
(S wave approximation). Such approximation is recovered in our formalism requiring £j = 
for all pairs. These results are quoted in the second column up to K = 10, being fully 
converged, and reproducing exactly the variational results 1 . The relevance of large relative 
orbital angular momenta can be judged by looking at the last column, where almost exact 
HH results up to K = 24 are given, the difference in some cases being larger than 200 MeV. 
This effect was not appreciated in Ref. [3] since their importance was estimated using only 
one gaussian for the radial wave function. 



Using the variational method of Ref. [37j, where nonzero relative orbital angular momenta 
were considered, we have obtained 3861.38 MeV for the energy and 0.363 fm for the root 
mean square radius, respectively, of the (L, S, I) = (0, 1, 0) ccnn four-quark state. This 
result is in perfect agreement with that obtained using the HH formalism: 3860.65 MeV and 
0.367 fm. 

We have also reproduced the calculation of the (L, S, I) = (0, 1, 0) ccnn state of Refs. 0, 



13| using the BCN model. In Table IIXI (second and third columns) we present the energies 



and RMS obtained for this state for all values of K. For K = 26 we have obtained an enere 



of 3899.2 MeV as compared to 3904.7 MeV of Ref. [13j and 3931.0 MeV of Ref. |9j. Ref. [13j 
designed a powerful method similar to the stochastic variational approach to study this 
particular system. Although they are not fully converged, the agreement gives confidence 
on both results. Ref. [9j uses a diagonalization in a restricted Hilbert space, obtaining a 
larger value. 

B. The QQnn system 

In Ref. the question: Does the quark model naturally predict the existence of QQnn 
bound states? was posed. The answer was clear, no compact bound four-quark state were 
found for any set of quantum numbers if only two-body potentials in a complete basis 



1 We have redone the calculation of Ref. [14] with the parameters used in this work for a proper comparison. 
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were used. One cannot discard that a modification of the Hilbert space, like for example 
considering only diquark configurations, or of the interacting hamiltonian, like many-body 
contributions, could give rise to bound states. Is the same conclusion still valid in the 
QQnn sector? To answer this question we have performed an exhaustive analysis of the 
QQnn spectra by means of th e q uark models described above. Some particular results 



were already presented in Ref. [3l|], where we studied the possibility of the X(3872) being 
a ccnn tetraquark. We gave arguments that favored the existence of QQnn stable states 
in nature, while making harder the existence of QQnn stable states. To make the physics 
clear we compared two particular sets of quantum numbers: J PC (I) = l ++ (0) for ccnn and 
J P (I) = l + (0) for ccnn. Based on a variational study with a confining mass independent 
many body potential we argued in Refs. [3l|, H(| that the binding would increase when 
increasing the mass ratio of flavor exotic four-quark systems. Definitive conclusions can 
only be obtained based on realistic calculations. For these purposes, in the present work 
we have considered all isoscalar and isovector states with total orbital angular momentum 
L < 1. For positive (negative) parity four-quark systems the ground state corresponds to 
L — (L — 1), since parity can be expressed in terms of the relative orbital angular momenta 
associated to the Jacobi coordinates as P = (— )* 1 +« ! +*3_ This means that P = —1 needs 
three units of relative orbital angular momentum to obtain L = [t\ = £2 = £3 = 1) while 
only one is needed for L = 1. The same reasoning applies for P = +1 states. Therefore, 
since the complexity of the calculation and the computing time increase with L, we have 
not considered L = 1 positive parity states that should be higher in energy. The calculation 
has been done up to the maximum value of K within our computational capabilities, K max . 

In Table IXl we present the ccnn CQC final results obtained for all possible L = isoscalar 
and isovector states and for the negative parity L = 1 states. We indicate for each state the 
maximum value of K used, K max . A first glance at this table sheds two general conclusions. 
Firstly, opposite to the QQnn case there exist bound states in the QQnn spectra when only 
two-body potentials are used. Second, once the four-body problem is properly solved the 
number of bound states is small. Curiously, there are two cases that do not converge to the 
lowest possible two-meson threshold but to a higher one, the (L = 0, / = 0) J p = 0~ and 
2~. These states get stuck in the lowest D wave threshold, either the D\ D 2 \d (5316 MeV) 
or the D*D*j,\ D (5325 MeV), satisfying A E > 0. 

As discussed in Sec. (Till the convergence of unbound or molecular states close to the two 
meson threshold is a difficult numerical problem. In order to minimize this problem we have 
introduced the extrapolation formula (|HJ). In Table IXIj we compare the energies obtained for 
K max and for K — > 00, using the extrapolation formula. For unbound or loosely bound states 
we observe how the extrapolation leads to four-quark energies within only a few MeV of the 
corresponding threshold. Throughout this manuscript we shall always refer to E iq (K maiy ) as 
the four-quark state energy, E iq . In particular cases we will use the extrapolation to study 
the characteristics of some specific states, mainly those that are close below the threshold. 

There are two ingredients that may alter the stability of the QQnn system: either the 
mass of the heavy quark or the interacting potential. One should wonder if the characteristic 
spectrum obtained in Table |X] would be greatly influenced by them. It was pointed out 
in the early 80 's that a QQQ'Q' four-quark state should be stable against dissociation 
into QQ' + QQ' if the ratio vtlq / vtlqi is large enough |1| . This was corroborated by chiral 
perturbation theory and lattice QCD studies of the bbnn system [8, 39]. Trying to disentangle 
if a proliferation of states is predicted when the mass of the heavy quark augments, we studied 
all ground states of the bottom sector using the same interacting potential as above. The 
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results are presented in Table Kill We observe that all bound states become deeper than in 
the charm sector and a few new states appear. Our results strengthen the conclusion that 
the larger the ratio of the quark masses the larger the binding energies. 

The second ingredient to be tested is the dependence of the results on the particular choice 
of the quark-quark interaction. To answer this question we have reanalyzed all the predicted 
bound states using the BCN model. The results are summarized in Tables IXIIII and IXIVI 
The existence of bound states is also evident when the BCN model is considered, although 
the characteristics of each state depend on the model considered. Of particular interest is 
the observation that the bottom sector presents, independently of the quark model, bound 
states with binding energies of the order of 100 MeV that should be possible to observe. 

In connection with the interacting potential used, it has been recently analyzed in Ref. [401 
the stability of QQnn and QQnn systems in a simple string model considering only a mul- 
tiquark confining interaction given by the minimum of a flip-flop or a butterfly potential in 
an attempt to discern whether confining interactions not factorizable as two-body poten- 
tials would influence the stability of four-quark states. The ground state of systems made 
of two quarks and two antiquarks of equal masses was found to be below the dissociation 
threshold. While for the cryptoexotic QQnn the binding decreases when increasing the mass 
ratio rriQ/m n , for the flavor exotic QQnn the effect of mass symmetry breaking is opposite. 
Although more realistic calculations are needed before establishing a definitive conclusion, 
the findings of Ref. |40| corroborate our results. 



C. Charm sector 

In the following we will study more closely those ccnn quantum numbers, J P {L, S, I), that 
may host a bound state (A# < in Tables IXl and [XTTT]) trying to unveil states that might 
be a consequence of model and computational approximations. Expected bound states are 
summarized in Table IXVl 



1. 1+(0,1,0) 

The possible existence of a four-quark state with these quantum numbers was predicted 
more than twenty years ago by Zouzou and collaborators Since then, several works have 
been devoted to study these particular quantum numbers by means of different methods 
and interactions, either in the charm or in the bottom sector 



Table [IX] we show the results obtained in this work with the BCN and CQC models. In 
both cases the state converges below threshold. For the CQC model the predicted binding 
energy is large, — 76 MeV, and Ar < 1, what would fit into the defined compact states. 
Opposite to that, the BCN model predicts a rather small binding, —7 MeV, and the RMS 
is larger than one but still not converged, although not increasing linearly. This state would 
naturally correspond to a molecule. 

Although this state would be stable against dissociation into two mesons, it may de- 
cay electromagnetically or weakly. The electromagnetic transition (QQnn) — > (Qn) (Qn) 7 
would be allowed if E iq is larger than the mass of two D mesons. The process is illustrated 
in the upper part of Figure [3J For the CQC model the four-quark energy is below the 
D D threshold and therefore the predicted state could only decay via a second order weak 
process into either two kaons (Cabibbo allowed) or two light mesons (Cabibbo suppressed). 
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The process is illustrated in the bottom part of Figure [3j Opposite to that, the results of 
the BCN model allow for an electromagnetic decay with the emission of a photon with an 
energy lower than 127 MeV in the ccnn rest frame. 

2. (1,2, 3)~ (1,2,1) and 1~(1,0,1) 

Bound states have been obtained either with the CQC or the BCN models or with both. 
In these detailed analysis of the possible thresholds is required. As previously 

noticed, the interacting models only consider central terms, L and S being proper quantum 
numbers. However, as illustrated is Sec. IIII Al the thresholds may be different in the coupled 
or uncoupled schemes, the former one being the relevant when trying to compare with 
experiment. We show in the upper part of Table IXVIl the results obtained compared to the 
threshold in the uncoupled and coupled schemes. The differences are noticeable, when the 
coupled scheme is used both interacting models give results above the corresponding lowest 
threshold, what discard these quantum numbers as promising candidates for being observed 
exp er iment ally. 

As already discussed in Sec. HVBt increasing the mass of the heavy quark will favor the 
binding. This is illustrated in Fig. H] comparing the four-quark energy of the 1~(1, 0, 1) state 
to its uncoupled threshold in the CQC model. As can be seen the CQC model predicts a 
tiny binding energy, but the extrapolation gives a very stable value around A^ « —6 MeV. 
Opposite to the case of unbound states, the RMS does not grow linearly with K. We observe 
the changes in A^ when increasing tuq, noting how A^ crosses zero for masses slightly below 
the charm quark mass. As a consequence, all the states discussed in this section might be 
candidates to be stable in the bottom sector due to the binding gained with the heavy quark 
mass. 

D. Bottom sector 

We repeat the same analysis as before for the bottom sector. Expected bound states are 
summarized in Table IXVl 

1. 1+(0,1,0) 

Opposite to the charm sector, in the bottom sector both quark models offer the same 
prediction, a compact deeply bound four-quark state. Within the CQC model the binding 
energy gets a value of —214 MeV while BCN predicts a somewhat smaller value of — 144 MeV. 
Concerning the possible decays, both of them are also below the threshold for electromagnetic 
decays into BB (10588 MeV for CQC and 10602 MeV for BCN) and therefore the only 
allowed decay mode for this state will involve a second order weak process into two open- 
charm (Cabibbo allowed) or two light mesons (Cabibbo suppressed). 

2. 0+(0,0,0) and 3~(1, 2,1) 

The existence of a positive parity QQnn bound state with quantum numbers S = and 
1 = was proposed in Ref. [B| using the BCN model if the ratio between the masses of the 
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heavy and light quarks was larger than 5—10, i.e., bcnn and heavier. We have obtained a 
similar limit for the CQC model, m&/m n ~ 15. 

For the J p = 3~ state a bound state has also been found using both models, being 
A^ = —140 MeV and A^ = —119 MeV, respectively. As can be seen in Table IXVIl it is the 
only member of a multiplet below all possible thresholds and therefore all strong decays are 
forbidden. 

No bound states were observed in the charm sector with these quantum numbers, being 
therefore consequence of the binding gained due to the larger heavy quark mass. Both 
states would present electromagnetic decays, the former one (bbnn) — > (B B\p) 7, with a 
photon energy in the range of ^ 400 MeV, and the latter through (bbnn) — > (B B\d) 7 or 
(bbnn) — > (B* B*\ s ) 7 emitting a photon of the order of < 400 — 450 MeV. 



3. 1-(1,0,0) 

Although no bound state is observed in the ccnn spectra, a bound state, A^ = —11 MeV, 
appears in the bottom sector with the CQC model. The structure of the wave function, 
« 1.182, clearly points to an extended meson-meson molecule instead of a compact 
four-quark state. Concerning its possible decay channels, being this state below all possible 
thresholds it could only undergo a weak decay. 



4. 1-(1, 0, 1), 2+(0, 2, 1) and (1, 2)~ (1, 2, 1) 

All these states are predicted to be bound in the uncoupled scheme but none of them 
survive the coupled thresholds. Although interesting from the theoretical point of view to 
carefully test our calculating framework, they are not expected to be observed in nature. 



E. Beyond QQnn 

Once the QQnn (Q = c, b) states have been discussed, few avenues remain to be explored. 
Among them, the consideration of QQss and QQ'nn states would be of interest. Concerning 
the strangeness ±2 states only those containing b quarks seem likely candidates to accom- 
modate bound states. We have redone the calculation for the most promising states, those 
with quantum numbers I(J P ) = 0(1 + ) and Q = c or b. In both cases the four-quark system 
is above the corresponding threshold. This can be easily understood because the mass ratio 
Mq I m q has diminished due to the large mass of the strange quark, thus increasing the contri- 
bution of the kinetic energy. Besides, the attractive one-pion exchange does not contribute. 
Our conclusions coincide with Ref. the binding energy decreases when decreasing the 



mass ratio Mg/nig. The second system, QQ'nn is made of distinguishable heavy quarks, 
i.e., bcnn among others, and therefore demands a modification of the current formalism, 
namely consider all possible combinations of £1, Su- This would lead to an increase in the 
basis size and consequently in computing time. However, one can draw a simple conclusion 
in light of Ref. [lf|, or the first six lines of Table III. The binding energy increases with the 
reduced mass of the heavy quark pair. Thus, one would expect a bound state with quantum 
numbers l + (0, 1,0) for the bcnn system, while the others appearing in the doubly-bottom 
sector would not be clear candidates for bound states with a be heavy diquark. 
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F. Meson- meson probabilities 



We present in Table IXVIII the meson-meson probabilities for some selected four-quark 
states according to Eqs. (fT3l) . These calculations where done by means of the variational 



method, Ref. [37|. As can be seen there is a perfect agreement with the HH results but the 
variational formalism allows to evaluate the probabilities of the different physical components 
in a simpler manner. Unbound states converge to two isolated mesons, the lowest threshold 
of the system, its RMS and Ar being very large. In contrast, bound states have a radius 
smaller than the threshold and they present probabilities different from zero for several 
physical states, the lowest two-meson threshold being contained in the physical four-quark 
system. Such states would be called compact in our notation. When the binding energy 
approaches the threshold, the probability of a single physical channel converges to one, what 
we defined as a molecular state. 



G. BCN vs. CQC 

The delicate interplay between the OGE and the chiral interactions in the description 
of the hadron spectra and baryon-baryon interaction has been widely discussed in the lit- 



erature 42]. In the QQnn spectra the relative strength of these interactions is even more 
important since A# is very sensitive to any modification in the chiral/OGE rate. The reason 
is that the threshold is not affected by the chiral part of the interaction since it is made of 
two Qn mesons, and there are no boson exchanges between heavy and light quarks. This is 
not the case for the QQnn systems, where bosons may be exchanged between the two light 
antiquarks. Therefore, an arbitrary modification in the strength of the chiral part of the 
interaction in any four-quark state where it is attractive, could bind the system. 

We therefore emphasize the importance of testing any model against as many observables 
as possible in order to constraint its parameters and ingredients. In this respect the heavy- 
light four-quark states are ideally suited for this task. Since only the total energy, and not 
the threshold, depends on the boson exchanges, the comparison between the predicted and 
measured Ae would provide us with precise information regarding the role played by these 
interactions in the meson spectra. 

Amazingly, as observed in Table IXVI almost all bound states are predicted independently 
of the dynamics. In other words, whenever a bound state is found with the CQC model, the 
BCN potential predicts a similar state. In general predicted binding energies are smaller for 
gluon-based interactions. Only the 1~(1, 0, 0) bbnn state is predicted by one of the models, 
the CQC. To illustrate the larger binding predicted by models considering boson exchanges 
we have selected one of the states whose binding energy in the uncoupled scheme has been 
found to be smaller than 10 MeV, the bbnn 2 + (0, 2, 1). Once boson exchanges are switched 
off the state lacks of enough attraction to be bound, behaving like two isolated mesons. 
These results are illustrated in Table IXVIIII To illustrate the difference in the structure 
induced by the boson exchanges we have plotted in Fig. [5] the evolution with K of Ar in 
both cases. In the unbound case the system is separating very rapidly while the bound one 
starts to converge to a large, but finite, value. 
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V. EXPERIMENTAL OBSERVATION 



The most promising mechanism for the production of a four-quark state is the indepen- 
dent formation of two cc and qq pairs to later on merge into a ccqq state. Following this 
idea the rate of ccqq production was estimated in Ref. [12j as R CC qq/ ' R ccq ~ 1/10. Therefore, 
any facility able to produce double charm baryons in sizable quantities should be able, in 
principle, to observe four-quark states. 



So far, double charmed baryons have only been observed by SELEX Collaboration [43 
although its limited statistics, ~ 50 events, makes doubtful that four-quark states can be 
produced in quantities large enough to be statistically significant. In Ref. 0| the production 
rate of double charmed baryons by the COMPASS experiment was estimated to be of the 
order of 10 4 to 1.7 • 10 4 events, this will indicate that up to 1500 four-quark events may 
be produced by COMPASS 200 GeV proton beam 45j. The number of ccq and bhq events 



at Tevatron was estimated in Ref. 461 ]. obtaining values of 10 5 and 10 4 events respectively. 
This would yield values of the order of 10 4 four-quark events to be produced in the second 
run of Tevatron. Although so far Belle Collaboration has not reported the observation of 



double charmed baryons, their production cross section was estimated in Ref. [47J, where 
values of the order of 10 4 events/year were obtained. This will translate into 1000 four- 
quark events/year. Theoretical predictions about the production cross-section of double 
charmed baryons at BaBar has been estimated, giving diverse values from hundreds to 



tens of thousands of events |48|. This will indicate that BaBar could produce up to 5000 
four-quark/events. On this respect it is worth to mention that theoretical cross-section 
predictions for double cc production have been found to be one order of magnitude too 



low as compared with experimental data |49j. This will imply that the predicted rate 



for double charmed baryon production may have been underestimated and therefore the 
expected number of four-quark events could be larger. In Ref. 0], the ccqq production rate 
at different facilities was estimated assuming that the dominant mechanism for double charm 
production at high energy colliders would be a disconnected double gluon-gluon fusion, 
(.9 + 9) + (9 + 9) ~^ ( c + c) + (c + c) . The predicted number of four-quark events produced at 
LHC by either LHCb or ALICE and at Tevatron is found to be very large, 9700, 20900, and 
600 events/hour respectively, while for RHIC they expect a smaller value of 12 events/hour. 

The picture that seems to emerge from these estimations is that nowadays we are on 
the verge where actual experimental facilities may be capable of start disentangling the 
properties of ccqq bound states. If this is not so, the future facilities that will be in operation 
in the next decade, if not sooner, at CERN and Fermilab will be able to provide a definitive 
answer to the existence of four-quark flavor-exotic states in nature. 



VI. SUMMARY 

In this work we have performed a systematic analysis of all ccnn and hhnn ground states 
within the framework of the hyperspherical harmonic method. In order to distinguish be- 
tween unbound, compact, and molecular four-quark states we have considered two different 
quark models widely used in the literature. We have analyzed both isoscalar and isovector 
systems with S=0, 1, and 2. We have considered L = states with positive and negative 
parity and L = 1 negative parity ones. The relevance of a careful analysis of the numerical 
thresholds together with the numerical approximations involved has been emphasized in 
order to avoid the misidentification of bound states. Estimations about the possibility to 
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detect these states on the next generation of experimental facilities have been performed. 

Our results are summarized in Table IXVI We have found five four-quark states that 
should be narrow and therefore possible to detect. Four of them are predicted independently 
of the interacting potential used, either CQC or BCN. The 1~(1,0,0) state is found only 
with the CQC model. The ccnn system would only have one bound state while up to four 
could be stable in the bbnn system. All predicted states are compact, only the 1~(1,0, 0) 
being molecular. Unfortunately, only one of them, the ccnn state with quantum numbers 
l + (0, 1, 0), is within the scope of the experimental facilities that will be available in the near 
future. 

Theoretical models point out to the existence of a double charmed isoscalar four-quark 
bound state with quantum numbers J p = 1 + , its properties depending on the quark model 
considered. The experimental detection and analysis of four-quark double charmed states 
will undoubtedly prove to be an invaluable testing ground to severely constraint the different 
theoretical models and therefore, it will allow to refine the theoretical predictions all over 
the hadronic spectra. 
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TABLE I: Energy difference, in MeV, between four-quark states QQnn {Q = c or b) and the 
corresponding two-meson threshold for different theoretical approaches. 
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a Results for two different set of parameters C\ (C2). 
^Results for two different interacting potentials. 
c Results for two different interacting potentials. 
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TABLE II: D and B meson masses (in MeV) and root mean square radius, RMS, (in fm) obtained 
with the quark models described in Sec. HH Experimental masses (Exp.) are taken from Ref. 291 ] . 
except for the state denoted by a dagger that has been taken from Ref. 30]. 
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"According to the Particle Data book [29( this signal can be interpreted as stemming from several narrow 
and broad resonances in the range 5650—5750 MeV. No quantum numbers are given. 
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TABLE III: Experimental lowest two-meson thresholds for charmed and bottom four-quark states 
in the uncoupled (UN) and coupled (CO) schemes as defined in the text. M\M 2 \l indicates the 
lowest threshold and E(M±, M 2 ) its energy. Energies are in MeV. 
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10650 


BB*\ P 


10604 


2"(1) 


D*D*\ D 


4317 


DD*\ P 


3877 


ID* D* 1 

d n \p> 


11023 


BB*\ P 


10604 


3"(0) 


D* D*\ P 


4018 


D* D*\p 


4018 


B*B*\ P 


10650 


B* B*\ P 


10650 


3"(1) 


D*D*\ D 


4317 




4290 


t> rs \D 


11023 


B B\\d 


10977 
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TABLE IV: Same as Table [TTT1 for ccnn states with the CQC model. We have evaluated the RMS, 
in fm. 





UN 


CO 


J p (I) 


M l M 2 \ L 


E(Mi,M 2 ) RMS 


M 1 M 2 \ L 


E(M 1 ,M 2 ) RMS 


L = 


0+(0) 


D l D\ P 


4426 


0.589 


D 1 D\ P 


4426 


0.589 


0+(l) 


DD\ S 


3872 


0.440 


DD\ S 


3872 


0.440 


1+(0,1) 


DD*\ S 


3937 


0.454 


DD*\ SiD 


3937 


0.454 


2+(0) 


D* D*j\p 


4499 


0.607 


DD*\ D 


3937 


0.454 


2+(l) 


D*D*\ S 


4002 


0.468 


DD\ D 


3872 


0.440 


0"(0) 


D 1 D^p 


4980 


0.738 


DD*\ P 


3937 


0.454 


o-(i) 


D}D}\p 


4996 


0.746 


DD*\ P 


3937 


0.454 


1"(0) 


Di D*j\ P 


4988 


0.742 


DD\ P 


3872 


0.440 




Di D}\p 


4988 


0.742 


DD*\ P 


3937 


0.454 


2- (0,1) 


D}D*j\p 


4996 


0.746 


DD*\ P 


3937 


0.454 






L 


= 1 S = 









l-(0) 


DD\ P 


3872 


0.440 


DD\ P 


3872 


0.440 


1-(1) 


DD x \ SiD 


4426 


0.589 


DD*\ P 


3937 


0.454 






L 


= 1 S = 


1 






o-(o,i) 


DD*\ P 


3937 


0.454 


DD*\ P 


3937 


0.454 


1"(0) 


DD*\ P 


3937 


0.454 


DD\ P 


3872 


0.440 


1-(1) 


DD*\ P 


3937 


0.454 


DD*\ P 


3937 


0.454 


2- (0,1) 


DD*\ P 


3937 


0.454 


DD*\ P 


3937 


0.454 






L 


= 15 = 


2 






1"(0) 


D* D*\ P 


4002 


0.468 


DD\ P 


3872 


0.440 


1-(1) 


D*D*j\ s ,d 


4499 


0.607 


DD*\ P 


3937 


0.454 


2~(0) 


D* D*\ P 


4002 


0.468 


DD*\ P 


3937 


0.454 


2-(l) 


D*D*j\ s ,d 


4499 


0.607 


DD*\ P 


3937 


0.454 


3~(0) 


D* D*\ P 


4002 


0.468 


D*D*\ P 


4002 


0.468 


3-(l) 


D*D}\s,D 


4499 


0.607 


DiD\ D 


4426 


0.589 
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TABLE V: Same as Table IIVI for bbnn states. 





UN 


CO 


J p {I) 


M X M 2 \ L 


E(M 1 ,M 2 ) RMS 


M 1 M 2 \ L 


E(M 1: M 2 ) RMS 


L = 


0+(0) 


B\ B\p 


11101 


0.373 


B\ B\p 


11101 


0.373 


0+(l) 


BB\ S 


10588 


0.284 


BB\ S 


10588 


0.284 


1+(0,1) 


BB*\ S 


10612 


0.287 


B B*\s,d 


10612 


0.287 


2+(0) 


B*B}\p 


11128 


0.377 


BB*\ D 


10612 


0.287 


2+(l) 


B*B*\ S 


10636 


0.291 


BB\ D 


10588 


0.284 


0"(0) 


B\ B\\p 


11614 


0.462 


BB*\ P 


10612 


0.287 


o-(i) 


B J B j\P 


11620 


0.464 


BB*\ P 


10612 


0.287 


l-(0) 


B\ Bj\p 


11617 


0.463 


BB\ P 


10588 


0.284 




Bi Bj\ P 


11617 


0.463 


BB*\ P 


10612 


0.287 


2- (0,1) 




11620 


0.464 


BB*\ P 


10612 


0.287 






L 


= 15 = 









l-(0) 


BB\ P 


10588 


0.284 


BB\ P 


10588 


0.284 


1-(1) 


B B\ \s y D 


11101 


0.373 


BB*\ P 


10612 


0.287 






L 


= 15 = 


1 






o-(o,i) 


BB*\ P 


10612 


0.287 


BB*\ P 


10612 


0.287 


i-(o) 


BB*\ P 


10612 


0.287 


BB\ P 


10588 


0.284 


1-(1) 


BB*\ P 


10612 


0.287 


BB*\ P 


10612 


0.287 


2- (0,1) 


BB*\ P 


10612 


0.287 


BB*\ P 


10612 


0.287 






L 


= 15 = 


2 






1"(0) 


B*B*\ P 


10636 


0.291 


BB\ P 


10588 


0.284 


1-(1) 


£> £>j\s,d 


11128 


0.377 


BB*\ P 


10612 


0.287 


2"(0) 


B*B*\ P 


10636 


0.291 


BB*\ P 


10612 


0.287 


2"(1) 


t?* t?* i 
£> £>j\s,d 


11128 


0.377 


BB*\ P 


10612 


0.287 


3"(0) 


B*B*\ P 


10636 


0.291 


B*B*\ P 


10636 


0.291 


3"(1) 


1 

t> £>j\S,D 


11128 


0.377 


B\ B\d 


11101 


0.373 
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TABLE VI: Same as 


Table HV 


for the BCN model. 




UN 


CO 


J p {I) 


M l M 2 \ L 


E{M l ,M 2 ) RMS 


M 1 M 2 \ L 


E(M 1 ,M 2 ) 


RMS 


L = 


0+(0) 


D X D\ P 


4341 


0.544 


D X D\ P 


4341 


0.544 


0+(l) 


DD\ S 


3772 


0.424 


DD\ S 


3772 


0.424 


1+(0,1) 


DD*\ S 


3906 


0.447 


DD*\ s ,d 


3906 


0.447 


2+(0) 


D*D*\p 


4511 


0.577 


DD*\ D 


3906 


0.447 


2+(l) 


D*D*\ S 


4040 


0.470 


DD\ D 


3772 


0.424 


0"(0) 


D 1 D^p 


4910 


0.664 


DD*\ P 


3906 


0.447 


o-(i) 


D*D*\ P 


4982 


0.684 


DD*\ P 


3906 


0.447 


l-(0) 


Di D*j\ P 


4946 


0.674 


DD\ P 


3772 


0.424 




Di D*j\ P 


4946 


0.674 


DD*\ P 


3906 


0.447 


2" (0,1) 


D*jD*j\ p 


4982 


0.684 


DD*\ P 


3906 


0.447 






L 


= 1 S = 









1"(0) 


DD\ P 


3772 


0.424 


DD\ P 


3772 


0.424 


1-(1) 


DD x \ SiD 


4341 


0.544 


DD*\ P 


3906 


0.447 






L 


= 1 S = 


1 






o-(o,i) 


DD*\ P 


3906 


0.447 


DD*\ P 


3906 


0.447 


l-(0) 


DD*\ P 


3906 


0.447 


DD\ P 


3772 


0.424 




DD*\ P 


3906 


0.447 


DD*\ P 


3906 


0.447 


2- (0,1) 


DD*\ P 


3906 


0.447 


DD*\ P 


3906 


0.447 






L 


= 15 = 


2 






l-(0) 


D* D*\ P 


4040 


0.470 


DD\ P 


3772 


0.424 


1-(1) 


D*D*j\ s ,d 


4511 


0.577 


DD*\ P 


3906 


0.447 


2-(0) 


D* D*\ P 


4040 


0.470 


DD*\ P 


3906 


0.447 


2"(1) 


D*D*j\ s ,d 


4511 


0.577 


DD*\ P 


3906 


0.447 


3-(0) 


D* D*\ P 


4040 


0.470 


D*D*\ P 


4040 


0.470 


3-(l) 


D*D*j\ s ,d 


4511 


0.577 




4341 


0.544 
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TABLE VII: Same as Table |V] for the BCN model. 





UN 


CO 


J p {I) 


M X M 2 \ L 


E(M 1 ,M 2 ) RMS 


M 1 M 2 \ L 


E(M 1: M 2 ) RMS 


L = 


0+(0) 


B\ B\p 


11113 


0.355 


B\ B\p 


11113 


0.355 


0+(l) 


BB\ S 


10602 


0.282 


BB\ S 


10602 


0.282 


1+(0,1) 


BB*\ S 


10651 


0.288 


B B*\s,d 


10651 


0.288 


2+(0) 


B*B}\p 


11176 


0.364 


BB*\ D 


10651 


0.288 


2+(l) 


B*B*\ S 


10700 


0.294 


BB\ D 


10602 


0.282 


0"(0) 


B\ B\\p 


11624 


0.428 


BB*\ P 


10651 


0.288 


o-(i) 


B J B j\P 


11652 


0.434 


BB*\ P 


10651 


0.288 


l-(0) 


B\ Bj\p 


11638 


0.431 


BB\ P 


10602 


0.282 




Bi Bj\ P 


11638 


0.431 


BB*\ P 


10651 


0.288 


2- (0,1) 




11652 


0.434 


BB*\ P 


10651 


0.288 






L 


= 15 = 









l-(0) 


BB\ P 


10602 


0.282 


BB\ P 


10602 


0.282 


1-(1) 


B B\ \s y D 


11113 


0.355 


BB*\ P 


10651 


0.288 






L 


= 15 = 


1 






o-(o,i) 


BB*\ P 


10651 


0.288 


BB*\ P 


10651 


0.288 


i-(o) 


BB*\ P 


10651 


0.288 


BB\ P 


10602 


0.282 


1-(1) 


BB*\ P 


10651 


0.288 


BB*\ P 


10651 


0.288 


2- (0,1) 


BB*\ P 


10651 


0.288 


BB*\ P 


10651 


0.288 






L 


= 15 = 


2 






1"(0) 


B*B*\ P 


10700 


0.294 


BB\ P 


10602 


0.282 


1-(1) 


t?* t?* i 
£> £>j\s,d 


11176 


0.364 


BB*\ P 


10651 


0.288 


2"(0) 


B*B*\ P 


10700 


0.294 


BB*\ P 


10651 


0.288 


2"(1) 


t?* t?* i 
£> £>j\s,d 


11176 


0.364 


BB*\ P 


10651 


0.288 


3"(0) 


B*B*\ P 


10700 


0.294 


B*B*\ P 


10700 


0.294 


3"(1) 


1 

t> £>j\S,D 


11176 


0.364 


B\ B\d 


11113 


0.355 



TABLE VIII: Comparison among different numerical approaches to the ccnn system. Energies are 
in MeV. 



(L,5,I) Ref. [14] HH(£i4 = 


0) HH 


(0,0,1) 4155 


4154 


3911 


(0,1,0) 3927 


3926 


3860 


(0,1,1) 4176 


4175 


3975 


(0,2,1) 4195 


4193 


4031 
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TABLE IX: BCN and CQC (L, S, I) = (0, 1, 0) ccnn results. 





BCN 


CQC 


K 


E 4q (MeV) RMS (fm) 


E 4q (MeV) RMS (fm) 





4100 




0.310 


4109 0.314 


2 


3999 




0.326 


3990 0.320 


4 


3954 




0.345 


3931 0.331 


6 


3933 




0.364 


3903 0.341 


8 


3921 




0.382 


3887 0.348 


10 


3914 




0.398 


3878 0.354 


12 


3910 




0.414 


3872 0.358 


14 


3907 




0.428 


3868 0.361 


16 


3904 




0.441 


3866 0.363 


18 


3903 




0.453 


3864 0.365 


20 


3901 




0.464 


3862 0.366 


22 


3900 




0.474 


3861 0.367 


24 


3900 




0.484 


3861 


26 


3899 




0.492 




DD*\ S 


3906 




0.447 


3937 0.454 


A E 




-7 




-76 


Ar 




> 1 




0.808 
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TABLE X: ccnn CQC energy, E 4q = E iq (K max ) (in MeV), RMS (in fm), A R , and A E (in MeV) 
as defined in Eqs. (|7|) and ([10]) . M\ M 2 \l indicates the lowest threshold and E(M\, M 2 ) its energy 
as obtained from Table HV1 



K 


J P (L,S,I) 


MiM 2 | L 


E{M 1} M 2 ) 




RMS 4(3 




A E 


28 


0+ (0,0,0) 




4426 


4441 


0.624 


>1 


+15 


28 


0+ (0,0,1) 


DD\ S 


3872 


3905 


0.752 


>1 


+33 


24 


1+(0,1,0) 


DD*\ S 


3937 


3861 


0.367 


0.808 


-76 


24 


1+(0,1,1) 


DD*\ S 


3937 


3972 


0.779 


>1 


+35 


30 


2+ (0,2,0) 


D* D*j\ P 


4499 


4526 


0.987 


>1 


+27 


30 


2+fO 2 1) 


D* D*\q 


4002 


4025 


0.879 


>1 


+22 


25 


0- (0,0,0) 


D x D 1 \ P 


4980 


5374 


0.738 


>1 


+394 


25 


- (0,0,1) 


D*jD*j\ p 


4996 


5012 


0.982 


>1 


+16 


25 


1- (0,1,0) 


Di D*j\ P 


4988 


5021 


0.944 


>1 


+33 


25 


1- (0,1,1) 


D x D*j\ P 


4988 


5018 


0.982 


>1 


+30 


25 


2 _ (0,2,0) 


D*jD*j\ p 


4996 


5387 


1.237 


>1 


+391 


25 


2- (0,2,1) 


D*jD*j\ p 


4996 


5025 


0.991 


>1 


+29 


21 


1- (1,0,0) 


DD\ P 


3872 


3938 


0.726 


>1 


+66 


23 


1-(1,0,1) 


D Dx\s t D 


4426 


4426 


0.527 


0.894 a 


+0 


21 


(0,1,2)"(1,1,0) 


DD*\ P 


3937 


3996 


0.739 


>1 


+59 


21 


(0,1, 2)- (1,1,1) 


DD*\ P 


3937 


4004 


0.814 


>1 


+67 


21 


(1,2,3)^(1,2,0) 


D*D*\ P 


4002 


4052 


0.817 


>1 


+50 


19 


(1,2,3)"(1,2,1) 


D* D}\s !D 


4499 


4461 


0.465 


0.766 


-38 



a In this case the radius has still yet to converge for K = K max . Its extrapolation gives a value larger than 
one but finite. 
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TABLE XI: ccnn CQC energies and A#, in MeV, evaluated for K mSLX and using the extrapolation 
([8]) in the limit K — > oo. 



rP It a t\ 

J {L,b,I) 


T? f V \ 
^iq{^raa,x) 


As 


E iq {K = CO) 


As 


0+ (0,0,0) 


4441 


+15 


4429 


+3 


0+(0,0,l) 


3905 


+33 


3862 


-10 


1+(0,1,0) 


3861 


-76 


3856 


-81 


1+ (0,1,1) 


3972 


+35 


3914 


-14 


2+ (0,2,0) 


4526 


+27 


4501 


+2 


2+ (0,2,1) 


4024 


+22 


3991 


-9 


_ (0,0,0) 


5374 


+394 


5323 


+343 


o-(o,o,i) 


5012 


+16 


4983 


-13 


I" (0,1,0) 


5021 


+33 


4993 


+5 


1-(0,1,1) 


5018 


+30 


4992 


+4 


2~ (0,2,0) 


5387 


+391 


5307 


+311 


2-(0,2,l) 


5025 


+29 


5000 


+4 


i-(i,o,o) 


3938 


+66 


3865 


-7 


1-(1,0,1) 


4426 


+0 


4420 


-6 


(0,1, 2)- (1,1,0) 


3996 


+59 


3927 


-10 


(0,1,2)-(1,1,1) 


4004 


+67 


3918 


-19 


(1,2,3)-(1,2,0) 


4052 


+50 


3993 


-9 


(1,2,3)-(1,2,1) 


4461 


-38 


4461 


-38 



TABLE XII: Same as Table IXl for the bbnn system. 



K 


J P (L,S,I) 


M 1 M 2 \ L 


E(M l ,M 2 ) 




RMS 4g 


Ar 


A E 


30 


0+(0,0,0) 


B B\\p 


11101 


10952 


0.328 


0.881 


-149 


26 


0+ (0,0,1) 


BB\ S 


10588 


10606 


0.365 


>1 


+ 18 


22 


1+(0,1,0) 


BB*\ S 


10612 


10398 


0.220 


0.765 


-214 


24 


1+(0,1,1) 


BB*\ S 


10612 


10623 


0.310 


>1 


+ 11 


28 


2+(0,2,0) 


B*jB*\ P 


11128 


11144 


0.627 


> 1 


+16 


30 


2+ (0,2,1) 


B*B*\ S 


10635 


10636 


0.314 


>1 


+1 


23 


1-(1,0,0) 


BB\ P 


10588 


10577 


0.335 


>1 


-11 


19 


i-(i,o,i) 


B B\ \s y D 


11101 


10980 


0.276 


0.740 


-121 


21 


(0,1, 2)- (1,1,0) 


BB*\ P 


10612 


10650 


0.493 


>1 


+38 


21 


(0,1, 2)- (1,1,1) 


BB*\ P 


10612 


10666 


0.517 


>1 


+54 


21 


(1,2, 3)- (1,2,0) 


B*B*\ P 


10635 


10677 


0.483 


>1 


+42 


19 


(1,2, 3)- (1,2,1) 


73* 73* 1 

B Bj\ s ,d 


11128 


10988 


0.276 


0.732 


-140 
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TABLE XIII: Same as Table |X] for the BCN model. In all cases the value obtained for the radius 
has still yet to converge for K = K max . 





M l M 2 \ L E(M u M 2 ) 


E 4q RMS 4q A R A e 


26 1+ (0,1,0) 

21 1-(1,0,1) 
21 (1,2, 3)- (1,2,1) 


DD*\ S 3906 
DDi\s, D 4341 
D*D*j\ s>D 4511 


3899 0.492 > 1 -7 
4380 0.640 > 1 +39 
4502 0.492 0.853 -9 



TABLE XIV: Same as Table Km] for the bbnn system. 



K 


J P (L,S,I) 


M 1 M 2 \ L 


E(M U M 2 ) 




RMS 4(? 


A R 


A E 


28 


0+(0,0,0) 


B B\\p 


11113 


11061 


0.334 


0.941 


-52 


28 


1+(0,1,0) 


BB*\ S 


10651 


10507 


0.220 


0.764 


-144 


28 


2+(0,2,l) 


B* B*\ s 


10700 


10723 


0.353 


>1 


+23 


19 


1- (1,0,0) 


BB\ P 


10602 


10639 


0.433 


>1 


+27 


19 


I" (1,0,1) 




11113 


11037 


0.265 


0.746 


-76 


19 


(1,2, 3)- (1,2,1) 


B Bj\ s ,d 


11176 


11057 


0.264 


0.727 


-119 



TABLE XV: Summary of bound states. 



Quark content J p 


(L,5,J) Model 


Decay mode 


ccnn 


1+ 


(0,1,0) 


CQC 


Weak 








BCN 


Electromagnetic 


bbnn 


1+ 


(0,1,0) 


CQC 


Weak 








BCN 


Weak 




3 - 


(1,2,1) 


CQC 


Electromagnetic 








BCN 


Electromagnetic 




0+ 


(0,0,0) 


CQC 


Electromagnetic 








BCN 


Electromagnetic 




1- 


(1,0,0) 


CQC 


Weak 



TABLE XVI: Ae, in MeV, in the uncoupled (A^ N ) and coupled schemes (A^,°) for different states 
with the CQC and BCN models. See text for details 





CQC 


BCN 


J p 


[L,S,I) 


MiM 2 | L5 


aUN 


M X M 2 \ J 


aCO 
^E 


M X M 2 \ J 


aUN 
E 


M X M 2 \ J 


aCO 
^E 


ccnn, l - 


(1,2,1) 


4499 


-38 


3937 


+524 


4511 


-9 


3906 


+596 


2- 


(1,2,1) 


4499 


-38 


3937 


+524 


4511 


-9 


3906 


+596 


3" 


(1,2,1) 


4499 


-38 


4426 


+35 


4511 


-9 


4341 


+161 


1" 


(1,0,1) 


4426 





3937 


+489 


4341 


+39 


3906 


+474 


bbnn 1 _ 


(1,2,1) 


11128 


-140 


10612 


+376 


11176 


-119 


10651 


+406 


2" 


(1,2,1) 


11128 


-140 


10612 


+376 


11176 


-119 


10651 


+406 


3" 


(1,2,1) 


11128 


-140 


11101 


-113 


11176 


-119 


11113 


-56 


1- 


(1,0,1) 


11101 


-121 


10612 


+368 


11113 


-76 


10651 


+386 
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TABLE XVII: Probability of the lowest threshold Pm 1 m 2 f° r different bound and unbound four- 
quark states. Energies in MeV and radii in fm. 



Quark content J P (L,S,I) 


E/t q 


A E 


RMS A R 


Pm 1 m 2 


ccnn 


0+(0,0,l) 


3877 


+5 


30.49 60.29 


1.00 




1+(0,1,0) 


3861 


-76 


0.37 0.81 


0.50 


bbnn 


0+(0,0,0) 


10948 


-153 


0.33 0.89 


0.25 




1+(0,1,0) 


10397 


-217 


0.22 0.77 


0.50 



TABLE XVIII: Energies and RMS of the bbnn 2+(0, 2, 1) as a function of K for the CQC model 
with (A) and without (B) boson exchange potentials. 





A 


B 


K 


E iq (MeV) RMS (fm) 


E 4g (MeV) RMS (fm) 





10763.0 


0.2131 


10793.3 


0.2173 


2 


10701.8 


0.2201 


10740.6 


0.2250 


4 


10670.6 


0.2283 


10712.2 


0.2337 


6 


10658.0 


0.2360 


10700.3 


0.2425 


8 


10650.1 


0.2436 


10692.2 


0.2521 


10 


10645.6 


0.2509 


10687.2 


0.2628 


12 


10642.6 


0.2580 


10683.5 


0.2743 


14 


10640.7 


0.2647 


10680.7 


0.2874 


16 


10639.4 


0.2712 


10678.6 


0.3025 


18 


10638.5 


0.2775 


10676.6 


0.3203 


20 


10637.8 


0.2836 


10674.9 


0.3415 


22 


10637.3 


0.2897 


10673.3 


0.3670 


24 


10636.8 


0.2957 


10671.7 


0.3973 


26 


10636.5 


0.3021 


10670.1 


0.4325 


28 


10636.2 


0.3077 


10668.4 


0.4720 


30 


10636.0 


0.3136 


10666.8 


0.5146 


B*B*\ S 


10635.5 


0.2906 


10635.5 


0.2906 
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FIG. 3: Possible electromagnetic (upper part) and weak (lower part) decays for ccnn bound states. 
The electromagnetic decay goes through an intermediate virtual meson (gray box) whose quantum 
numbers will depend on those of the initial four-quark state, virtual D* for the case of J p = 1 + . 
The weak decay showed corresponds to the Cabibbo allowed one into two kaons that incorporates 
two vertex (c + u) — > (s + d) or (c + d) — > (s + u). Other diagrams involving two W emission into 
fermion-antifermion, (ccqq) — > {sq) + (sq) + (ff) + (ff), or mesons, (ccqq) — > (sq) + (sq) + (ud) + (ud), 
may also contribute [41]. 
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FIG. 4: Ae as a function of the heavy quark mass, mq, for 1 - (1, 0, 1) with the CQC model. The 
grey bands take into account the differences between the values obtained for K max , band upper 
part, and the result obtained with the extrapolation, black squares. 
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FIG. 5: Ar as a function of K for the 2 + (0, 2, 1) bbnn with boson exchange potentials (solid line) 
and without (dashed line) for the CQC model. 




32 



